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Abstract 

The theory of passive scalar transport in two dimensional turbulent fluids is 

generalized to the case of 2D MHD. Invariance of the cross correlation of scalar 

concentration and magnetic potential produces a novel contribution to the concentration 

flux. This pinch effect is proportional to the mean potential gradient, and is shown to 

drastically reduce transport of the passive scalar across the mean magnetic field when 
RmBq/ ^7ip{y 2 ^ > 1. Transport parallel to the mean magnetic field is unchanged. 

Implications for models of transport in turbulent magnetofluids are discussed. 
PAC NOS. 47.25.Jn, 47.65.+a 
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Recent work on the theory of magnetic potential (flux) diffusion in two 

dimensions 1 and the a-effect driven dynamo in three dimensions 2 has focused on the 

effects of small scale magnetic fluctuations. The thrust of the research pursued by several 

groups is that Lorentz forces produced by rapidly amplified small scale magnetic fields 

strongly "back-react" on the fluid dynamics, thus inducing a marked departure of mean 

field evolution from kinematic predictions. Thus, in two dimensions the magnetic potential 

diffusivity is determined by the competition between small scale hydrodynamic and 
magnetic energies, 3 i.e. r\ e g ~ ^v 2 ) - ^B 2 ) • Similarly, the a -effect is set by the imbalance 

of fluid helicity and magnetic hyper-helicity, 4 i.e. a ~ (w - a) - (b-Tj . The magnitudes of 

small scale magnetic energy and hyper-helicity are constrained by topological conservation 

laws which govern magnetic field evolution. Specifically, in two dimensions, conservation 

of mean square magnetic potential directly ties the magnetic fluctuation energy to the net 

spatial transport of mean magnetic flux. Similarly, in three dimensions, conservation of 
magnetic helicity directly relates the magnetic fluctuation hyper-helicity (^B ■ to the net 

mean-field a-effect. In both cases, the mean-field transport coefficient in question (i.e. the 
flux diffusivity in 2D, a-coefficient in 3D) is reduced by a factor (l+R m BQ/ / 4^:p^v 2 )) 

relative to the analogous kinematic prediction. Hereafter, the mean magnetic field-induced 
suppression factor (l+7? m .B 2 ) / / 4;zp^v 2 )) is written as (l+R) -1 , where 

R = R m B%l ^Tzp^v 2 ) • The notation is standard, with R m the magnetic Reynolds number. 

It follows that inhibition of magnetic flux diffusion and saturation of the dynamo process 
occur for mean field magnetic energies drastically below the equipartition level (i.e., for 
bI/4k> p(v 2 ) /R m )- Note that for R > 1, D e ff = Dk/(1 + R) scales proportional to the 

collisional resistivity, even if the turbulence correlation time is short relative to resistive 
diffusion times, i.e. % c r\l?-«\. This reflects the "freezing in" property of magnetic 
potential dynamics. Such results call for a re-evaluation of the conventional wisdom 
developed using kinematic mean-field electrodynamics. 5 In particular, the impact of the 
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novel mechanism for saturation of the dynamo at modest field levels upon predictions 
derived from conventional dynamo theory should be assessed. 

From the discussion above, one is naturally lead to wonder about the effects of 
small scale magnetic fields on turbulent transport and mixing in MHD fluids, in general. 
This question is obviously relevant to understanding the origin and profile of the solar 
differential rotation and the structure of the solar convection zone as well as to other 
problems in astrophysical fluid dynamics. Note that recent analysis of helioseismological 
data suggests that the solar convection zone is enmeshed by a network of small scale 
magnetic flux tubes with significant volumetric packing fraction, 6 thus suggesting that the 
convection zone should be viewed as an MHD fluid. The dynamics of each are governed 
by turbulent transport of momentum and heat, respectively. As small scale magnetic 
fluctuations are known to modify the effective (eddy) viscosity coefficients, correlation 
times and spectra of MHD turbulence, it is indeed reasonable to explore the effect of 
magnetic turbulence on turbulent transport in magneto fluids. Also, understanding the 
dynamo in a Reversed Field Pinch plasma requires a theory of the turbulent transport of 
magnetic helicity. 7 As turbulent momentum and heat transport in MHD fluids are 
exceedingly complex phenomena, we first seek to identify and understand the dynamics of 
a simple prototypical model. The most obvious such "hydrogen atom problem" is that of 
transport of a passive scalar in 2D magnetohydrodynamics (MHD), which is a 
straightforward generalization of the classic paradigm of passive scalar convection by a 
turbulent fluid. 8 This problem is contained within the simple, familiar model consisting of 
the magnetic potential, vorticity and passive scalar concentration advection equations in two 
dimensions, i.e. 



dy/ 



+ V0x£- Vi/a = rjV 2 y/, 



(la) 



V 2 + V0x£- VV 2 0- W 2 V 2 = vyx£- vvV, 



(lb) 
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— c + V0Xz-Vc-DoV 2 c = O. (lc) 



<p and i/A are the velocity stream function and magnetic potential, respectively. Thus fluid 

velocity v = V0 x £ and magnetic field B = Vy/ x £ . Here the collisional resistivity is 77 , 

the kinematic viscosity is v and the collisional diffusivity of the scalar concentration is Do • 

Hereafter, we take 77 = V for simplicity, but allow Dq^t\. This system has the familiar 

quadratic conserved (up to diffusive dissipation) quantities, which are total energy 
Jj 2 Jc|(V0) 2 + (V)/A) 2 jy2, mean square magnetic helicity \d 2 x{lff 2 )^ an d cross-helicity 

\ d 2 x(V(/) ■ Vy/) . It also has two other quadratic inviscid invariants, namely \d 2 x{c 2 l'l) 

and the (important) concentration-magnetic potential (C-MP) correlation \d 2 xcy/ . Note 

that the conservation of C-MP correlation is a consequence of the fact that both c and yi 

are conserved along fluid element trajectories. Thus c = c(0) , yf = y/{<p) so c = c(ys). 

For a given turbulent velocity field, a standard approach to the determination of the 

passive scalar flux F is to proceed via a quasi-linear closure. Specifically, we calculate the 

mean field flux of scalar concentration as 



r = ((V0x£)c) (2a) 



and approximate T as 



E = {( V<p c W) + (V 0W c)} x z (2b) 

where and 0^ are determined from Eqns. (lc) and (lb) respectively, assuming finite 
correlation time for the MHD turbulence. Thus, taking (c) and (iff) to refer to the mean 

concentration profile and magnetic flux (where (B) = V(i/a) x £), it follows that the 
responses c$ and 0^ are: 
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4 1) = - T ^(Y^x£.V(c), 



(3a) 



V 2 4 1) = T cfe Y(^xfv(vV W ) +(Y^ 1) x£-VV 2 (va) 



(3b) 



Contributions from the vorticity advection nonlinearity vanish, as {<pc) k = . Neglecting 
terms of order l/k 2 Ly/, where Ly is the characteristic scale length of the mean magnetic 
potential, one then finds: 



4 1) = T cfe(YY^x£-V(vA) 



(3c) 



Taking the 0, y/ and c fields to be isotropic then yields: 



k 



{W) k Y{c)-(Yc-Y¥) k Y{¥) 



(4a) 



= -D K V(c) + VV(y/) 



(4b) 



Note that the total scalar concentration flux consists of the usual diffusive piece -D^V(c), 
(with the usual kinematic diffusion coefficient Dk = £\Y0 / T ck) an d a novel magnetic 



potential gradient driven flux VV(y/) (with pinch velocity coefficient 

V = Y,(Yc'Yf}i ( T ck)- Since the V(i/a) flux vanishes along contours of constant mean 

k 

magnetic potential, only the cross-field flux is renormalized by the pinch term, i.e. 
En = E|| K , the kinematic result. 

The readily evident proportionality of the pinch V to the correlation (Vc- Vi/a), 

which is obviously related to the C-MP correlation, gives a clue as to the importance of the 

(topological) conservation of mean-square magnetic potential for passive scalar transport. 
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The conservation of C-MP correlation in the presence of mean profiles (c) and ((//) 
implies: 



4 (c W) + (W) • Y(c) + (yc) ■ V(va) = -{rj + D )( Yc ■ V» . (5a) 
at 



Noting that by definition (vi/a) = -rj^V^), (vc}j_ = -Dj^jf V(c) and requiring 

stationarity, the following relation between the evolution of C-MP correlation, transport of 
magnetic potential and scalar concentration, and the collisional dissipation of these 
quantities follows directly 

(%ff + D±_ eff )v(c) ■ V(V> = fa + £>o)(Yc ■ Vva) . (5b) 

Here 77 = D^/[l + ^], and recall D# is the kinematic turbulent diffusion coefficient. 
Thus, the correlation (Vc ■ Vi/a) may be obtained from Eqn. (5b) and substituted into Eqn. 
(4) for V , which yields 



fa + A)) 



Now, we are concerned with cross-field transport only, so we take V(c) parallel to V(i/a). 
Thus 

(v(y)) 2 / \ 

r± = - D K V± (c) + x c [D^ff + Veff) Vi (c) . (7) 

Use of the definition r_i_ = ~D± e ff Vj_(c) and some straightforward manipulation then 
finally reveal the effective cross field diffusivity as predicted by mean field theory to be 
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(l+R) 



1+^^R 
ri+Do 

l + ^-R 
. rj+Do 



(8) 



The form of D± e ff obtained above is non-trivial. For 77 = Do which corresponds 
to Prandtl number (as well as magnetic Prandtl number) unity, D±_ e ff = Dk(^+R)~ 1 - 
Hence, the cross-field diffusivity is suppressed, in comparison to the kinematic prediction, 
by the same factor as the magnetic potential diffusivity is reduced. This may be 
understood by noting that in the limit of V = 77 = Do , c - c(</>) and y/ = so c = c(y/) 

on all inertial scales. As magnetic potential transport is suppressed for 
R = ByRm/ p(v 2 ) > 1 > an d the scalar concentration is frozen into the magnetic potential, it 

is not surprising that the cross-field scalar diffusivity is similarly quenched. This 
interpretation is supported by the observation that for Do » 7?, D± e ff —> Dk, i- e - the 

kinematic result is recovered. The reversion to simple kinematics in this limit follows from 
the breakdown of the relation c = c(y/) due to the disparity between Do and 77 . Simply 

put, since the passive scalar concentration c is not frozen into the fluid, the fact that 
magnetic potential is frozen into the fluid (thus suppressing magnetic diffusion) is 
dynamically irrelevant. Curiously, for Do — > 0, D±_ e ff = DkI (1+^) 2 • In this limit, cross 

field passive scalar transport is quenched (for R> 1) even more strongly than magnetic 
potential diffusion is. Care should be taken in drawing conclusions from this result, 
however, since Do — > may not be compatible with a fully stationary state. Finally, the 
absence of D\\ renormalization may be understood by realizing that the suppression of 
magnetic flux transport is due to a kind of "elastic memory" of magnetic field lines in a high 
magnetic Reynolds number fluid. Thus, a "plucked" field line tends to "snap back" to its 
original configuration. Hence, the restoring force implied by the elastic memory analogy 
acts in a direction perpendicular to the mean magnetic field. Since c = c(y/), only the cross 

field transport is quenched. 
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The predictions presented above can easily be tested by a simple numerical 
experiment in which a mean magnetic potential profile (iff(x)) and a circular blob of 
concentration density c(x) centered at the origin are initialized in a turbulent 2D 
magnetofluid. The magnetic potential amplitude should be adjustable, so that cases with 

i?<l and R>1 can be compared. For R<1, the concentration blob will spread 
isotropically with (8r 2 )~ Dr? ■ For R>1, the blob should develop a pronounced 
anisotropy with variances along the principal axes given by \Sy 2 j~DK^ an d 
(Sx 2 ) ~ D±_ e ffT ~ Dkt/(1 + R) (for 77 ~ Do)- Keep in mind that for R m » 1 (the case of 
interest) (B) 2 /An «p(y 2 ^/2, so that the mean field can be well below the equipartition 

level. This is because the dynamically relevant field is the small scale magnetic field, the 
energy of which is amplified by stretching of (B). Here, the parameter R should be 
thought of as a measure of the strength of small scale magnetic field effects in terms of 
mean field strength. 

The least-defensible approximation made in this paper is that of a single correlation 
time x c k for each of the spectra ((f 2 ^ , (v /2 ) fc an d { c2 ) k - ^ n strong MHD turbulence, 

when inertial range equipartition occurs, it is however reasonable to expect both (v 2 ) fc and 

(b 2 )^ to have the same correlation time, which will, however, depend on Brms and Bq 

(i.e. via the Alfven effect). 9 Since c is a passive scalar, with dynamics determined by 
(\ 2 ) k and (indirectly) by (b 2 )^' it does indeed appear reasonable to hypothesize a single 

correlation time for this system. Another assumption made ab-initio here is that the 
microscopic collision frequency v c of the gas exceeds the (electron) cyclotron frequency 
(v c >Q e ), so that the passive scalar diffusion tensor is invariant with respect to the 
direction of Bq. If Q. e >v c , the collisional flux of certain scalar quantities becomes 

anisotropic. Specifically, electron heat conduction along the field has diffusivity 
X\\ = Vl/v e ,e but X L = V ee P 2 e - Moreover, DqV 2 T -> V M ^T M V M 2" + X^T, where 
V|| = (i?o + B ■ Vj j\B\. Hence, anisotropic transport introduces new non-linearities and a 

different transport mechanism mediated by magnetic fluctuations, only. 10 
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The implications of these results must be discussed keeping in mind the obvious 
question of the relevance of a 2D MHD model (where is conserved) to 3D dynamics 

where magnetic helicity, not mean square potential, is the relevant topological invariant. 

First, it should be noted that magnetized incompressible MHD fluids are well described the 
reduced MHD model, 11 where (V 2 ) invariance is broken only by the bending of B z (i.e. 

dl^ff^jdt ~ Bzlydijf/dzj. Thus, situations of approximate "mean" square potential 

conservation can arise particularly when r c k z V a < 1 (i- e - turbulence correlation time is 
shorter than the parallel Alfven transit time). The transport suppression predictions given 
here are directly relevant to such situations. Moreover, the theory discussed here could be 
extended to incorporate the effects of Alfvenic radiation in the mean square magnetic 
potential budget. Second, many aspects of the phenomena discussed here are generic to 
MHD turbulence. For example, mean flow evolution by turbulent transport is governed by 
the competition between fluid and magnetic stresses i.e. 

d ^ = d 
dt dxj 

Hence, amplification of ^ to equipartition levels will certainly affect (and likely reduce) 

turbulent momentum transport. This scenario, which closely resembles that of the 
prototypical model discussed above, is not encompassed by existing theories of solar 
differential rotation. Indeed, the latter treat only the effects of mean magnetic fields. 12 
Finally, the results suggest a re-evaluation of the lore concerning turbulent transport in 
magnetofluids. Specifically, the effects of small scale magnetic fields on eddy conductivity 
and mixing length theory should be investigated. Also, greater attention should be given to 
mechanisms for the spontaneous amplification of magnetic fluctuations, such as the 
magnetic shearing instability. 13 Here, quasilinear predictions of saturation levels with 



V;V. 



Anp ( 



(9) 



9 



(b 2 ^ l%7i > p (v 2 )/2 are intrinsic to the instability dynamics, thus avoiding the 

cancellations induced by growth to equipartition. 
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